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Mail from Albert Rich: 

Hello Josef, 

From your recent email to Michel, I was glad to see that the DNL is alive and well. You must 
be near celebrating its 25th anniversary. 

Since you have expressed an interest in Rubi in the past, I thought you might want to check 
out her new look at http://www.apmaths.uwo.ca/~arich/ 

Click on the "Integration Rules" menu item to find and view the appropriate set of rules that 
apply to a particular integrand type. For example, try finding the trig rules for integrands of 
the form (a+b tan(c+d x)) A n. Eventually I hope to continue this cascading navigation menu 
tree right down to the specific rule needed to integrate a given expression; thereby mirroring 
on the website the decision tree Rubi uses to integrate expressions... 

Please let me know what you think of the site and any suggestions may have. 

Aloha, 

Albert 


Leam more about “ Rubi ” m the next DNL, Josef. 

Mail from Francisco Marcelo Fernandez: 

Dear Derivians, 

I had a problem calculating a simple integral with Derive. It 
is explained in the attached file. May be I am missing something but 
I cannot find what. 

Best regards, Marcelo 


#1: [a :e Real (Q s ™) s m :e Integer (0, «)] 

Derive fails to calculate 


■XI 

f 2 ■ m 2 

#2: J x ■ EXP(- a- x ) dx 

—00 


#3: 


a? 


e 

— o a 


z 

a- x Z - m 

- x dx 


But it can do it with a little help 


00 

f Z ■ m Z 

#4: Z-J x ■ EXP(- a-x ) dx 

0 


- (2-m + l)/2 ( l'] 

V " tJ ! 


#5: 
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Dear DUG Members, 

It is end of June now and I am very 
happy that I can offer DNL#98 in time. 

The “Truth Tables ..." is one of the 
many “Contributions waiting to be pub- 
lished" from page 2 which found its way 
from my stack of more or less ready 
papers into this newsletter. MacDonald 
Phillips sent this article long ago (2007) 
and I am very sorry that I didn't pre- 
sent it earlier. I hope that Don will ac- 
cept my excuse. 2007 was the time 
when TI-92, Voyage 200 and TT-89 
were full in use while TI-NspireCAS was 
in its infancy. DERIVE does not need a 
special function or program because 
TRUTH_TABLE() does a great job. TT- 
NspireCAS is another case: there is no 
function implemented and we cannot 
create own drop down menus. So it was a 
challenge for me to provide a truth ta- 
ble tool for TT-NspireCAS. Pages 30 
and 31 show the result. 

It has happened very often that ar- 
ticles come in which are of immediate 
interest (see Michel's note on the solu- 
tion of cubics in this issue) or that I 
become fascinated by a special problem 
(mail, journal, talk, ...) which then gets 
its self dynamics. 


can, started DERIVE, and then I was 
caught by programming, browsing my 
books, searching the web, ... and the 
result is a 10 pages article. I hope that 
you will enjoy the treatment of the 
"Misunderstood Sister of the Fibonacci 
Sequence". 

Michel Beaudin sent a toolbox for 
TI-NspireCAS wich has its history in its 
DERIVE past. Michel had provided DE- 
RIVE - tools for his students at ETSMTL 
several years ago. As he and his stu- 
dents changed to TT-Nspire he updated 
this useful collection for the new tech- 
nology and added some functions which 
are implemented in DERIVE. All func- 
tions, programs and written instructions 
and provided examples are in French. I 
asked Michel for permission to produce 
an English version which was given with- 
out hesitating. Many thanks, Michel. 

I will meet Michel and a couple of 
other DUG-members in July in Kalamata, 
Greece, at ACA 2015. 1 am sure that we 
will have an interesting meeting despite 
the economic problems which Greece 
and its people must face in these days. 

I will give a report in DNL#99. 
Always Yours 


Take as an example the Stern- 
Brocot Sequence treatment in DNL#98. 
I read about it in the Scientific Ameri- 


Josef 


Download all DAX-Derive- and Tl-files 



http : / / www . austromath . at/ dug/ 


In December we will celebrate DNL#100 (= 25 Years DUG). It will 
be great if some of you - especially members from the early DUG 
Years - will contribute for this very exceptional issue. All articles, 
notes, comments, memories, ... are very welcome, Josef. 


EDITORIAL 


DNL 98 


P 2 


I 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE , 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 


Next issue: September 2015 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encryption, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Liike-Rosendahl, GER 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

A New Approach to Taylor Series, D. Oertel, GER 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 

Rational Hooks, J. Lechner, AUT 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 
Statistics of Shuffling Cards, H. Ludwig, GER 
Charge in a Magnetic Field, H. Ludwig, GER 

and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Solving Cubics 

Michel Beaudin, Montreal, Canada 

Hello my friends : if you have time, please take a look at the PDF « Lost_something ». It is about the 
solutions to the following third degree polynomial equations: x 3 - 3x + 1 = 0, x 3 + 3x + 1 = 0, 
x 3 + 4x 2 + x - 1 and 27x 3 + 216x 2 + 584x + 528 = 0. I have used the following CAS : Derive, Maple, 
Nspire CAS and WolframAlpha Pro. 

The first one and third one have three real distinct roots while the second one and fourth one have one 
real root and two complex conjugates. 

In the case of 3 distinct real roots, Derive (version 6.1, the last one), is using trig expressions for the 
answers and this is very clean and clear. And when Derive is using Cardano’s formula, NO radicals 
appear in the denominator. 

After the Derive screen images, you will find the same four examples done with Maple 15 (I don’t 
have the latest version). Then screen images from Nspire CAS (OS 4, latest version) will show that 
Nspire CAS can do a better job with some help : I have defined some functions in my TI-Nspire CAS 
library « kit ETS MB » and the results are similar to those of Derive. I am proud of this. 

Solutions by WolframAlpha were printed and scanned and are shown on the pages following the solu- 
tions by Nspire CAS. I find easier to use WolframAlpha (Pro) instead of Mathematica... but I have 
Mathematica 7 on my computer because I need it to run Rubi! 

Best regards from Montreal, Canada, 

Michel Beaudin 


Derive 


#1: SOLVECx - 3-x + 1 = 0, x) 

#2: x = 2-COSl I v x = - 2-COSl I v x = 2-SINl 


{-T-) ” ’ ■ 2C0 S (t) v ‘ ■ 2SIB br) 


#3: SOLVECx + 3-x + 1 = 0, x) 

1/3 1/3 

(4-75 - 4) (4-75 + 4) 

#4: x = + L' 


1/3 


4 

1/3 


1/3 1/3 'l 

73.(4.75-4) 73- (4.75 + 4) 

+ 


(4.75 - 4) (4.75 + 4) 

+ + L- 


1/3 1/3 'i 

73.(4.75-4) 73. (4.75 + 4) 

+ 


1/3 1/3 

(4.75 - 4) (4.75 + 4) 
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#5: SOLVE (x +4-x + x - 1 = 0, x) 

5.^3 


#6: 2 • 7l3 • COS 

X = 




2-713-SIN 


’ » 


3 J 


2-713 -SIN 


4 ™("“T") 1 


#7: 


#8: 


4 
3 

1 2 
SOLUTIONS(27 • x + 216-x + 584-x + 528, x) 

1/3 1/3 1/3 

2.(789 - 9) 2- (789 +9) 8 (789 - 9) 


1/3 


9 3 

1/3 'l 


73.(789 - 9) 73.(789 + 9) 

+ 


(789 - 9) 


1/3 


1/3 

(789 +9) 8 


1/3 

(789 +9) 8 


9 


1/3 1/3 ^ 

73.(789 - 9) 73.(789 + 9) 

+ 


9 


9 


Maple 


> solved — 3 -x + 1 = 0, x) ; 


|(-4 + 4,/T) ,,3 + - 2-—,-i ( -4 + 4./T) 

(-4 + 4I/T) 


1/3 


(-4 + 4I/T) 


1/3 


1/3 


+ 1^3 (4- (-4 + 4IV3 ) 


1/3 


(-4 + 4I/T) 
1 

(-4 + 4I/T) 
(-4 + 4I/T) 173 


,-l(-4+4l/J) 


1/3 


YTj-yWs (4(-4 + 4I>/3) 


1/3 


> solved + 3 -x + 1 = 0, x) ; 


t ( 4+4 ^) 1/3 + t — ^ t ( 4+4 ^> 1/3 

( 4 + 4 /r) 


(4 + 4VT) 1/3 ' 2 ‘ V " l 2 

,4(4 + 4/V) 1/3 - 


Li/yf-±(4 + 4/y) 


1/3 


(4+4VT) 1 
-ll^f-l(4+4/T)'' 3 - 


(4+4V5) 


1/3 


(4 + 4/T) 


1/3 
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solve{x 3 + 4-x 2 +jc-1=0,jc); 


j ( ~ 260 + 156IVT) 173 + — 26 1/3 - J, 

0 3 ( -260 + 156 1 V 3 ) 

--jV (-260 + 156lVT) 1/3 - — — -J7J 

3 (-260 + 156lV 3 ) 

+ 4 1/3" [ 4 ("260 + 156lVT) 1/3 

2 6 


26 


3 (-260 + 156I/T) 
13 

3 (-260 + 156IVT) 

+ 156l/T) 1/3 - 


1/3 


1/3 


, — —■ ( -260 + 156 iVT) 


1/3 


4 - 4-iVT ( 4 (-260 


26 


3 (-260 + 1 56 1 V 3 ) 
> solve( 27*x A 3 + 216*x A 2 + 584*x + 528 = 0,x); 


(9 + /89-) 173 - 


1/3 


1 


(9 + V89 ) 


1/3 


9 (9 + V~89") 1/3 3 ’ 9 

2 -v+ivrf4(9+v^) 1/3 


9( 9 + /89-) 1/3 3 ' ‘ v ~ [ 9 

,-j{9+JW) 


9{9 + JW) 

2 

9(9 + V~89") 
2 

9 (9 + 7^9") 


1/3 


1/3 


1/3 


1/3 


-- 7 -^ 4 ( 9 + 789 ) 


1/3 


TI-N spire CAS and my contribution : equations x 3 - 3x + 1 = 0 and x 3 + 3x + 1 = 0. 
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TI-N spire CAS and my contribution : equation x 3 + 4x 2 + x - 1 = 0. 



1 1 3 _ , 2 , \\ "Error: Resource exhaustion 

exact lcZeros(27- v +216. i +584- .v+528,.v]j 

k 

it_ets_mb\compact ciibic\21 \ 3 +216- v“'+584- .r+528,.v) 
11 11 

| 2- (J89+9) 3 2- (J89-9) 3 8 ~(j89+9) 3 ( (^- 9 ) 3 8 ! 

i 1 

1 

. -(Js 9 +9) 3 * 

[9 9 3’ 9 9 3 

9 9 ( 

9 


& Wolfram Vlpha 




1 


H ■ « 1 ■ 

solvE(x"3-3xtl=fl ! Jt) 

B) 





|npj|- inlerpretation: 

salve x 


■R^ulLS: 




i+j vr 


i +fVT 


x - -- 


i vr i 

. — . i 

2 ^ 




i 

r > 

IT 




x 
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Input inter piretjt on. 

solve + 3 x + 1 * 0 for jt 


Results: 



input [nterpiefat^n: 

solve jt 3 * 4 jc 2 + x - 1 = 0 or * 


Results; 


^ « *3.65109 
jr*s -0.726109 



Roots in thf complex plane: 

tmtx) 

I.MlO' * 

O *■ R*00 

-I.KlO -10 
-2^I0’ IS 

V 

-I -2 -\ 0 
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Input interpretation: 

solve 27* 3 +2lBje 2 + 584* + 528 = 0 


Results- 


2 

jc = - 

9 


-- 


8 

3 


= -j+ Z(1 ~ ,V5 l -^(i-nV8)^9tVm 
3 aVstVw 9 

8 2 f 1 + f V 3 ) 1 / / — \ / jf ~~ 

= -r + , =--[i-w)p+-4& 

3 9^9 + 789 " 9 


Reactions on Michel’s paper: 

Albert Rich: 

Hello Michel, 

Yes, what is simpler may be in the eyes of the beholder, but Maple and Mathe- 
matica returning real solutions involving the imaginary unit certainly isn’t... 

Aloha from Hawaii, 

Albert 

David Jeffrey: 

Hello Michel, 

The trig formulae are usually called Viete's formulae. If a cubic has 3 real roots, 
then Viete's formulae give them without needing sqrt(-l). 

Congratulations on your simplification of Nspire formulae. I assume that your sim- 
plifications could be used on any trig expression, whether it came from a cubic or 
elsewhere. 

Radicals appear in the denominator of Cardano formulae, because in the general 
case this is necessary to prevent branch cuts jumping to an incorrect solution. In 
all cases it is clear that Maple simply substitutes in a general Cardano formula, 
rather than selecting the formula according to the case. 


David 
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David Stoutemyer: 

Vietes trig formula CAN involve arcsin or arccos of a quantity having abs greater 
than 1. See Wikipedia for a hyperbolic alternative for this case. 


Jose Luis Galan: 

Michel, congratulations for the simplification in Nspire. And, as in many situations, 
our old good friend Derive keep being the most appropriate !!! 

Best, 

Jose Luis 

Josef Bohm: 

many thanks to Michel for his impressive cubics' collection which I can accomplish 
with the WxMaxima and Mathcad 15 solutions. 

Thanks also for the additional comments of the CAS-experts. 

I'd like to put this in the next DNL's User Forum. 

I am trying to program the exact solution of quartics for Nspire but I must admit 
that I failed until now because of the many nested roots which cannot be handled 
by the Nspire - at least in my opinion - but I will proceed with my experiments. 

Best regards to you all, 

Josef 


(%il) 

<%ol) 


WxMaxima 14.12.1 


solve ( [x A 3-3*x+l=0 ] , [x] ] ; 



4 $ %i l ) 

2 2 } 




4/3 

+ 


'x/s"' 1 


2 2 




7 (%i2) solve ( [x A 3+3*x+l=0 ] , [x] > ; 


(%o2) 


Vs" %i 1 



-Vi 1 %i 1 


2 2 
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(%i3) solve ( [x A 3+4*x''2+x-l=0] , [x] ) ; 


13 ! 


rvr 


(%o3) [x- 


%i l'j 

2) 


{ 


' 13 %i 

G5\ 

.2 3 3/2 

54 J 


f 13 %i 65 N | 

n 

( t/T %i l'j 

Li 

f 13 %i 65'' 

Us 3/2 54 J 

l 2 2 ) 

3 7 x 

,2 3 3/2 54 , 


1/3 


f -Vi 1 Sri 1 I 

i'-jv %i _ 1 j + 13 ^ 2 2 ' 


9 1 ■ 


3 3/2 54 


4 ( 13 %i es'i 173 

x= — — 1 +- 


f 13 Sri 65V /3 3 ^2 3 3/2 54 J 


13 


V 


(%o4) [x= — 


[^4 


r a -vii n + fl '| 1/3 

v. 72 5 aij 


L 72 9 81 


rv 


’ 13 %i 

65V 

2 3 3/2 

54 j 

f 



8 

2 ) 

3 / * 


4 

--] 


81 


(L/T %i l'j 

L 8 

f Vi 1 Sri 1 

1 . , 

[ B VSS 1 , 8 "j 

l 2 2 

l 2 2 ) 

81 ( 

a Vis 1 _ a Y 

r Jf — . 

1/3 3 ! 

V 72 9 8 ij 

725 aij 


aVi^ a 1 1/3 3 

81 1 + — 

72 5 B1 


8 8 ) 
72 9 8 1/ 


] 


1/3 


Mathcad 15 


x - 3x + 1 auflosen 


- 1 


J.+&Y + 

2 2 ) 


1 4±i\ 

2 2 J 


_L_ 

3 

2 

3 


1 >/3'i^ r- f 1 V3*i^ . /— . 

- + 1 \ — yj 3’ — + I • i + ■y 3‘ i 


2 2 2 


2 2 j 


2 - 


1 , Vj-L 


h 


2 2 J 


r 1 V 3 -i V r r f 1 V 3 -iV 

— + -* — ! + 1 -/ 3 -i + / 3 - — + v 1 ■ 


2 2 j 


2 2 J 


•1 


1 4±±\ 

2 2 ) 


2 ■ 
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3 2 

27x + 216x + 584x+ 528 auflosen — » 


128-V89 10880 V 

59049 53144lj 


8-^89 _V| 

729 81 ) 


729 81 J 


I 3 

3 


.,h'12»jS.222»') -432 ♦ »V 5 * sihT^-^ + — 'I 

59049 53144 y ^ 729 81 j ^ 59049 531441 / 


729 81 J 


■.„„f 8 -v /89 8 ^ ( 128-^89 10880 ^ „ r- . r- ( 128-^89 10880^1 

432 - + — | + 81 - + , — 8 + 8 i-v 3 + 81 i-v 3 ’ + , 

^ 729 81 J \ 59049 531441 / ^ 59049 53144 V 


, 62 . 1*42 4 - 1 '] 

729 81 J 


Wilfried Zappe, Germany: Problems with normpdf on TI-Nspire 



*Dok^ 


0.75 


fl(x) =norm Pdf(x 



* 


- 0.63 .. 


Analysing the graph we can detect the maximum 
but not the inflection point. 


It is not possible to plot the first derivative of the 
distribution function. 




It seems to be that normpdf internally is not de- 
fined as a function. 

One has to define his/her own distribution func- 
tion. Then it works as you can see on the next 
screen shots. 
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Hubert Langlotz, Germany: Just another bug? 

While correcting the papers of our final examination I found another bug in OS 4.0. 

Its real bad, because the use of the tool “bounded area” was a possible way to solve this problem in the 
exam! 






DNL 

Same happens in version 3.9. Interestingly enough it works properly when introducing a 
slider for m. Then everything is ok. 
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The Mathematician and the Clockmaker 

Josef Bohm, Wiirmla, Austria 

In the German issue of Scientific American, Spektrum der Wissenschaft , from May 2015 I found an 
interesting article about a sequence: “The Misunderstood Sister of the Fibonacci Sequence” written by 
Jean-Paul Delahaye. 

This sequence is defined by s 0 =0, s x = 1, s 2n = s n , s 2n+l =s n + s n+v 

It is called Stern [1] -Brocot [2] Sequence. This sequence offers many very remarkable properties which 
are worth to be investigated. 

The following short program generates n elements of the sequence: 

st_bro(n P seq, el p k) := 

Prog 

seq := [1] 
k := 2 
Loop 

If k > n exit 

el := IF(M0D(k, 2) = Q P seqx(k/2) p seqi((k - l)/2) + seqxCfk - l)/2 + 1)) 
seq := APPEND(seq, [el]) 
k :+ 1 
seq 

st_bro(20) = [1, 1, 2, 1, 3, 2, 3, 1, 4, 3, 5, 2, 5, 3, 4, 1, 5, 4, 7, 3] 

There is no special structure visible. Things will change when we produce a graph presenting the ele- 
ments versus their position numbers in the sequence: 

st_bro_pts (n , k, seq P el p pts) := 

Prog 

seq := [1] 
pts := [ [1 P 1]] 
k := 2 
Loop 

If k > n exit 

el := IF(M0D(k s 2) = 0, seqj,(k/2) P seqiCCk - l)/2) + seqi((k - l)/2 + 1)) 
seq := APPEND(seq, [el]) 
pts := APPENDCpts, [[k, p el]]) 
k :+ 1 

pts 

st_bro_pts (10000) 

DERIVE makes it possible to plot the first 10 000 points: 


500 



) 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 9000 9500 10000 11 


-50 
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Zooming in we can find a fractal similar structure: 



6000 6200 6400 6600 6800 7000 7200 7400 7600 7800 8000 8200 8400 8600 


I will follow the article in the journal and form groups of the elements of size 2 k - omitting the 0. 
I store one sequence (5000 elements) as st_seq. 


st_seq := st_bro(5000) 

VECT0R(|" st_5eq ], k, 0, 5) 

L [2 A k 2 A (k + 1) - 1] J 


[1] 


. [ 1 , 6 , 5 , 





[1, 

2] 




[1, 

3, 

2, 

3] 

[1, 

4, 

3, 

5, 

2, 

5, 

[1, 5, 4, 7, 3, 

8, 

5, 

7, 

2, 

7, 

9, 4, 11, 7, 10, 3, 11, 8, 13, 5, 

12, 

7, 

9, 

2, 

9, 


3, 4] 

5, 8, 3, 7, 4, 5] 

7, 12, 5, IB, 8, 11, 3, 10, 7, 11, 4, 9, 5, 6] 


I don’t need the centred order of the groups of elements. So I produce the first 12 groups and select 
them one after the other with maximum 32 elements: 


grps := VECTOR(rst_seq ], k, 0, 12) 

L [Z A k Z A (k + 1) - 1]J 


[1] 

[ 1 , 2 ] 

[1, 3, 2, 3] 


[1, 

4, 

3, 

5, 2, 5, 3, 4] 








[1, 

5, 

4, 

7, 3, 8, 5, 7, 2, 7, 

5, 

8, 3, 7, 4, 5] 






[1, 

6, 

5, 

9, 4, 11, 7, 10, 3, 

11, 

8, 13, 5, 12, 7, 9, 

2, 9, 7, 12, 5, 13, 8, 11, 3, 

10, 

7, 

11, 

4, 9, 5, 6] 

[1, 

7, 

6, 

11, 5, 14, 9, 13, 4, 

15 

, 11, IS, 7, 17, 10, 

13, 3, 14, 11, 19, 8, 21, 13, 

IS, 

5, 

17, 

12, 19, 7, 16, 5 


[1, 8, 7 , 13, 6, 17, 11, 16, 5, 19, 14, 23, 9, 22, 13, 17, 4, 19, 15, 26, 11, 29, IS, 25, 7, 24, 17, 27, 10, 23 
[1, 9, 8, 15, 7, 20, 13, 19, 6, 23, 17, 28, 11, 27, 16, 21, 5, 24, 19, 33, 14, 37, 23, 32, 9, 31, 22, 35, 13, 3 
[1, 10, 9, 17, 8, 23, 15, 22, 7, 27, 20, 33, 13, 32, 19, 25, 6, 29, 23, 40, 17, 45, 28, 39, 11, 38, 27, 43, 16, 
[1, 11, 10, 19, 9, 26, 17, 25, 8, 31, 23, 38, 15, 37, 22, 29, 7, 34, 27, 47, 20, 53, 33, 46, 13, 45, 32, 51, 15 
[1, 12, 11, 21, 10, 29, 19, 28, 9, 35, 26, 43, 17, 42, 25, 33, 8, 39, 31, 54, 23, 61, 38, 53, 15, 52, 37, 59, 2 

The last list is created by: 


grps 


12 , 1 , [ 1 


32 ] 
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By the way, appending a 1 to all groups of numbers we create palindromes. 


We need one more step to find the property which I want to show. Look at the numbers of the first 
column, the second column, the third column, and so on: 


[1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 1] 


[2, 

3, 

4, 

5, 

6, 

7, 

8, 

9, 

10, 

11, 

12] 

[2, 

3, 

4, 

5, 

6, 

7, 

8, 

9, 

10, 

11] 


[3, 

5, 

7, 

9, 

11, 

13, 

15, 

17, 

19, 

21] 

[2, 

3, 

4, 

5, 

6, 

7, 

8, 

9, 

10] 




[5, S, 11, 14, 17, 20, 23, 26, 29] 

[3, 5, 7, 9, 11, 13, 15, 17, 19] 

[4, 7, 10, 13, 16, 19, 22, 25, 28] 

[2, 3, 4, 5, 6, 7, 8, 9] 

[7, 11, 15, 19, 23, 27, 31, 35] 

[5, 8, 11, 14, 17, 20, 23, 26] 

[8, 13, 18, 23, 28, 33, 38, 43] 

This seems not to be too spectacular: all columns form arithmetic sequences. Now notice the differ- 
ences of these sequences: 

[0, 1, 1,2, 1, 3, 2, 3, 1, 4, 3, 5, ...] - Ooops: 
st_bro(12) = [1, 1, 2, 1, 3, 2, 3, 1, 4, 3, 5, 2] 
this is again the Stern-Brocot Sequence. 

What’s about the sums of the groups? (You may remember that the sums in the rows of the Triangle of 
Pascal are the powers of 2.) 

VECTOR (Kg rps ), n, 12) 
ri,l 

[1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, 177147] 

FACT0R([1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, 177147], Number) 

[ 2 3 4 5 6 7 8 9 10 111 

Ll, 3, 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 , 3 ,3 J 

You may try to proof this result! (According to the author of the article it should not be too difficult.) 

VECTOR (MAXCgrps ), n, 12) 
n,l 

[1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233] 

I am sure that you can recognize these numbers! Yes, you are right, it is the Fibonacci Sequence. 

The above mentioned Triangle of Pascal is a very well known structure of numbers. We can discover 
two nice properties which are possibly not so well known. 
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pas Cp) := VECTORCEVECTORCCOMBCn, k), k, 0, n)L n, 0, p) 

It’s no problem to present the Triangle: 

We observe the numbers appearing in the 
increasing diagonals - the numbers of the pas (7) = 
same colours. 



ds := d lags (20, 20) 


[1] 

[ 1 , 1 ] 

[ 1 , 2 , 1 ] 

[1, 3, 3, 1] 

[1, 4, 6, 4, 1] 

[1, 5, 10, 10, 5, 1] 

[1, 6, IS, 20, IS, 6, 1] 

. [1, 7, 21, 35, 35, 21, 7, 1] . 


di ags (p , m) produces the first m increasing diagonal sequences of a Pascal triangle with p rows. 
These are the last 5 diagonals: 


[ 1 , 14 , 78 , 220 , 330 , 252 , 84 , 8 ] 

[ 1 , 15 , 91 , 286 , 495 , 462 , 210 , 36 , 1 ] 

[ 1 , 16 , 105 , 364 , 715 , 792 , 462 , 120 , 9 ] 

[ 1 , 17 , 120 , 455 , 1001 , 1287 , 924 , 330 , 45 , 1 ] 

. [ 1 , 18 , 136 , 560 , 1365 , 2002 , 1716 , 792 , 165 , 10 ] . 

Now let’s calculate the sum of the elements in the 20 rows of the above matrix (which are the diago- 
nals): 

VECTORCICds ), k, DIM(ds)) 
k,l 


[ 1 , 1 , 2 , 3 , 5 , 8 , 13 , 21 , 34 , 55 , 89 , 144 , 233 , 377 , 610 , 987 , 1597 , 2584 , 4181 , 6765 ] 
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We meet again the ... Fibonacci Numbers. Let’s count the odd numbers in the diagonals. You may 
start with the first coloured diagonals above. I ask again the CAS for doing the job: 

VECTOR(X(MQD(ds p 2)), k, DIM(ds)) 
k,l 

[1, 1- 1, 3, 2, 3 , 1, 4, 3, 5, 2 , 5, 3, 4, 1, 5, 4, 7 , 3] 


Do you recognize this sequence? Isn’t this funny? 

Can you imagine that there is a relationship between the binary representation of a number, continued 
fractions and the elements of the Stem-Brocot Sequence? I’ll show you. 

Let’s take the binary representation of 100 [3] : 

DUAL (100) = [1, 1, 0, 0, 1, 0, 0] 


I count the ones and the zeros how they appear in the list: 2, 2, 1, 2, and close with a 0 (for no one is 
following). Then I use the list [2, 2, 1, 2, 0] to form a continued fraction: 


0 + 


1 


1 

2 + 

1 

1 + 

1 

2 + 

2 


7 


19 


You might ask: “Ok, but where is the connection with the Stem-Brocot Sequence?” Wait just one 
moment. We look for elements s m and s m of the sequence. 

(st_bro(101)) = [7, 19] 

[100 p 101] 

xV ioo is the numerator and its successor s m the denominator of the fraction. To go for sure we’ll take a 
second example: 


DUAL (5 03 ) = [l p 1, 1, 1, 1, 0, 1, 1, 1] 

(st_bro(504)) = [23, 6] 

[503, 504] 


I wrote a small function for calculating the continued fraction and apply it immediately: 

fr(x, f, i) := 

Prog 

f := 1/xjT 
i := 2 
Loop 

If i = DIM(x) exit 
f := l/(f + x|i) 
i :+ 1 
xj.i + f 


23 

fr([5, I, 3]) = 

6 


I am quite sure that you didn’t expect this relationship. 
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We didn’t mention the most remarkable property of this sequence until now. It was the reason for 
Stem to introduce his sequence in 1858. 

We all know the usual way to count the rational numbers using the (1, 2, 3, . . .) x (1 ? 2, 3, . . .) table 
and counting the fractions in the diagonals. There are rational numbers which appear several times by 
cancellation (1/3 = 2/6 = 3/9 = . . .). Omitting these fractions disturbs the system and makes it difficult 
to find the fraction on say position n. 

The Stern-Brocot Sequence helps solving this problem. We calculate a new sequence 

*0 *1 *2 Vl 



ra t_n Limbs (n) := VECTOR 


(st_bro(n + 1)) 

k 


(st_bro(n + 1)) 

k + 1 


k p n 


rat_mjmbs(2Q) = 



3 

2 



4 

3 


3 

5 


5 

2 


2 

5 




3 

S 


This gives all rational numbers. And the n rational number in this sequence is the fraction - JL - L which 

S n 

can easily be calculated using the binary forms described above. Which is the 1000 th rational number 
in this list? 

DUAL (1000) = [1 P 1, 1, 1, 1, 0, 1, 0 P 0 P 0] 


11 

f r([5 p 1, 1, 3, 0]) = 

39 

It should be 1 1 ! Let’s check! 


(st_bro(1000)) = 11 

1000 


Remembering the palindromes formed by the 2 k - element groups of the original sequence it is not so 
astonishing that the sequence of fractions - all rational numbers in a nice order - will also appear in 
pretty groups with 2 k elements each: 


r_n := rat_ruimbs(200) 
VECTOR (f r_n 

L [2 A k 


2 A (k + 1) - 1] 


j, k, 0, 4) 
[ 1 ] 

1 


2 


L 2 


13 2 

, , , 3 

3 2 3 




The fractions form a binary tree which is called Calting-Wilf-Tree. 
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We can sort the elements in increasing order: 

VECTOR C|"S0RTCr_n )], k, 0, 4) 

L [2 A k Z A (k + 1) - 1] J 

[ 1 ] 

1 

, 2 

2 

2 

3 

3 

4 

4 

5 




8 7 


3 
2 

4 

3 

5 

4 ’ 5 


5 5 

! 

3 2 

7 8 


7 7 

4 ’ 3 


Compare with the graph of this binary tree (Reference [4]) 


1 

0 






A A 

1 2 3 3 

A 

4 5 

A 

5 4 

4 5 5 4 

AA AA 

1233 4554 
5 7 57 7 S T 5 

3 3 

A A 

5 7 57 
4 5 54 

2 I 

A A 

7S 7 5 

3 3 2 1 


8 1 

Performing a binary search say for — and starting at - , we have to move right, left, right and left 

again on the branches (= RLRL). Each finite string consisting of Rs and Ls defines exactly one ra- 
tional number. Irrational numbers are defined by infinite sequences of Rs and Ls. 


Task for the reader: 

Find out, which number is defined by the sequence RLRLRLRLRL .... 

And here is another one: RL° RLR 2 LRL 4 RLR 6 LRL 8 . . . (R 2 = RR and L 4 = LLLL etc.) 
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This binary tree can be constructed by calculating the mediant of two consecutive elements of the se- 

p p' p + p' 0 1 

quence. The mediant of two fractions — and — is . We start with fractions — and — . We 

1 f i f 1 A 

q q q+q 1 0 

know very well that the latter is no fraction but in this case it is helpful to take it as one. 

The mediant has a nice representation as the sum of two vectors: 



Let’s follow the “mediant - idea” with DERIVE : 

NUMERATQR(a) + NUMERATOR (b) 

medfa, b) := 

DENOMINATORS) + DENOMINATORS ) 



meds(frs, nfrs) 

Prog 

frs := REVERSE (REST (REVERSE (f rs ) ) ) 
nfrs := [FIRST(frs)] 

Loop 

If DIM(frs) = 1 exit 

nfrs := APPENDfnfrs, [medffrsxl, frsjZ), frsj.2]) 
frs := REST(frs) 

nfrs := APPENDfnf rs , [l/IF(nfrs = [0] p I, nfrsj.2), r, l/0 r ']) 


meds([Q, 1/0]) = [0, 1, 1/0] 


meds([Q, 1, 1/0]) = 


0 , p 1 , 2 , 1/0 

2 



1 


112 3 

meds 

0 , , 1 , 2 P 1/0 

= 

0 , , p , 1 , , Z P 3 P 1/0 


2 


3 2 3 2 


This is a good occasion to apply the ITERATES-command: 

ITERATES (meds (1), 1, [0, 1, 1/0], 2) 


[ 0 , 1 , 1 / 0 ], 


0 , , 1 , 2 , 1/0 

2 


112 3 

0 , , , , 1 , , 2 , 3 , 1/0 

3 2 3 2 


]] 


But we would like to see the binary tree growing. So let’s build a matrix of the vectors created above 
including the very first row: 


meds_rows(n) := APPEND([ [ [0 , l/0]]] p VECrGR([v] ? v, ITERATES (meds (1 ), l f [0, 1, 1/0] p n))) 
meds_rows (3 ) 
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[Q, 1/0] 

[0 P 1, 1/0] 

1 

0 , , 1 , 2 , 1/0 

2 


112 3 

0, , , , 1, , 2, 3, 1/0 

3 2 3 2 


1 

1 

2 

1 

3 

2 

3 

4 

3 

5 5 

Op 

, 

, 

, 

, 

, 

, Ip 

, 

p 2, 

4 

3 

5 

2 

5 

3 

4 

3 

2 

3 2 


It looks good, but it is not really the tree we would like to see, because the new rows are containing the 
elements of the previous rows, too. They need to be eliminated. I remember a DERIVE - command 
REMOVE_ELEMENTS - and this should do the job: 


bin_tree(n) := APPEND([[[Q, 1/0]]], VECTOR (TREM0VE_ELEMENT5((meds_rows(ri)) , (meds_rows(n)) )], k, n + 1)) 

L k + 1,1 k , 1 J 


|bin_tre&(4)| 



And this is really the binary tree from above which can be used for the binary search. 


I will close with another property of this really remarkable sequence: 

P 1 

The simplicity of a fraction — is defined as . I omit the first row of the bi n_tree generated 

q p-q 

matrix and calculate the simplicities of all elements: 




r 1 i 


VECTOR 

VECTOR 

, k, DIM(w) 

, w , VECTORS , v, REST(bin_tree(n))) 



DENOMINATORS ) ■ NLJMERATOR(w ) 

1 



l k k j 

/ 


simpls(4) 


[ 1 ] 


1 

2 


1 

2 


1111 
3 6 6 3 






' 1 

1 

1 

1 

1 

1 

1 

1 ' 







. 4 

10 

15 

12 

12 

15 

10 

4 . 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

5 

14 

24 

21 

28 

40 

35 

20 

20 

35 

40 

28 

21 

24 


1 

14 


1 


5 


Do you recognize a special feature of the matrix? Add the elements in the rows! 
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simp1s_sums(n) := VECTOR 


* 

r f 

I 

VECTOR 




DENOMINATORS ) ■ NUMERATORS ) 
k k 


k, DIM(w) 


, vi, VECTOR(v » v, REST(bin_tree(n))) 
1 


simpls_sjms(10) = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1] 

Finally Jean-Paul Delahay presents another very simple formula to generate an element r n+ \ of the 
sequence of all rational numbers from r n \ 

1 


" +1 1 + 2-Lr.J-r/ 

which gives us another opportunity to generate this sequence applying DERIVE ’ s ITERATES: 


St(x) := 


1 + 2 ■ FLOOR(x) - x 




f 

l 1 + 2 ■ FLGOR(x) - 


st_i ts(n) := ITERATES| » x, 1, n 

3 1_ its (20) 

1 132 1 4 3 5 2 5 3 1 5 4 7 3 8 

1 * i 2 , P p p 3 p p p p p p p P 4 p p p p p P 

2 5 23 43525 3 4 547385 


[1] Moritz Abraham Stem(1807 - 1894) was a German mathematician. He introduced this sequence 
1858 by inventing this way to count the rational numbers. 

[2] Achille Brocot (1817-1878) discovered this sequence - completely independent of Stern - by 
investigating transmission ratios of gearwheels. He was a clockmaker (see reference [4]). 

[3] This is a tool from an earlier DERIVE Newsletter 


I will close now presenting the Stern-Brocot Sequence and some of its properties. It was interesting - 
and sometimes a little bit challenging programming tools for this paper. Please note that there are a lot 
of resources in the web dealing with this “Misunderstood Sister of the Fibonacci Sequence”. 
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Truth Tables on the TI-89 
by MacDonald Phillips 
don.phillips@gmail.com 


Truth Tables on the TI-89 is a simple program to 
use. Simply put in the logical expression and 
press ENTER. First, however, mn the menu() 
program that is located in the Tmth folder. Under 
FI: Tmth Tables, is the main program and logical 
connectors. 

Besides the logical connectors included on the TI- 
89 (not, and, or, xor), I created two functions, one 
for implication and the other for if-and-only-if. 
IMP() takes two arguments, IMP(p,q), meaning if 
p then q. IFF() also takes two arguments, 
IFF(p,q), meaning if and only if p then q. 

Under F2: Vars is a listing of the common letters 
used as variables in symbolic logic and the last 
logical expression used. 


dHRBBfFTT 

,]K!: 


3 

noT 

and 

4 

or 

5 

xor 

6 

IMPC 

7 

IFFC 


TYPE DR USE + [ENTER] UR [ESC] 


[Truth Tabl«l 

jfJMf 'i 




# :p 

3sq 

4: r 
5: s 
6:t 
7: u 
8: o 


i 

TYPE UR USE + [ENTER] UR [ESC] 


The tmth table for implication is created thus: 


f f1t 

' F£t' 


Truth Tables 

Yars 



Truih-T<INP<p,q>> 

TRUTH RflD HUTU FUNC 0/30 


And the output is: 


You are given a list of the logical variables and 
the expression entered (exp) and the possible 
tmth values of the variables and resulting tmth 
value of the entered expression. Press ENTER to 
display the rest of the tmth table and/or to exit the 
program and return to the HOME screen. 


r fit 

|Truth Tables 

' F£t' 

Yars 

] 

[p q 

exp] 




[t t 

t] 




[t f 

f] 




[f t 

t] 




[f f 

t] 




TRUTH 


RhD HUTU 

FUNC 

— isron 


r fit 

|Truth Tables 

' F£t' 

Yars 

] 

jZ 

-i— 1 

■ 

tCinpCp , 

q» 

Done 

Trut,h_T 

:imi 

= <p>q 



I TRUTH 

RhD HUTU 

FUNC 

i/30 
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Some simple logical expressions are (unfortu- 
nately) automatically reduced to true or false by 
the TI-89. If that is the case, then either Tautol- 
ogy or Self-Contradiction will be displayed in- 
stead of the truth table. 


( f1t 

' F£t‘ 


Truth Tabl« 

Yars 



Truth_Ttp or not, p) 

truth 

RHD HUTD 

FUNC 

0/30 


( f1t 

' F£t‘ 



1 

|Truth Tables 

Yars 



1 


Tautology 


TRUTH RHP HUTD FUHl liillHt 


Complex arguments may also be evaluated by the 
program. For instance, evaluate the argument 
if(if(p or q then r) and not r then not p and not q). 


f f1t 

' F£t‘ 

1 

|Truth Tables 

Yars 

1 


This can be entered as Truth_T(IMP(IMP(p or q, 
r) and not r, not p and not q). 


Truth-T <IMP<IMP<P or g,r>... 

TRUTH RHP RUTD TUNC 0/30 


The truth table is: 


r fit 

|Truth Tables 

' F£t‘ 

Yars 

: ] 

[p q r 

exp] 

[t t t 

t] 


[t t f 

t] 


[t f t 

t] 


[iff 

t] 


[fit 

t] 


TRUTH 

RHD HUTD FUNC IHItH 


and 


As can be seen, the truth value of the argu- 
ment is true no matter what the true values of 
the variables; the argument is therefore a tau- 
tology. 


( f1t 
T ruth Tables 

' F£t‘ 

Yars 


[p q r 

exp] 

[f t f 

t] 


[f f t 

t] 


[f f f 

t] 


TRUTH 

RHD HUTD FUNC llilW 


(Note: The program will support any number of logical variables, but only about 7 or 8 can be dis- 
played in a row in the output screen. To display the truth table on the HOME screen just type in tt and 
press ENTER.) 

Reference: Logic: An Introduction , Lionel Ruby, J. B. Lippincott Company, 1960. 

...it is unworthy of excellent men to lose hours 
like slaves in the labor of computation. 

Gottfried Wilhelm Leibniz 
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It is much easier for DERIVE users because the function TRUTH_TABLE is ready made 
implemented. The logic operators AND, OR, XOR, NOT, IMP, and IFF (logic equivalence) 
are available. 


Let’s create the truth table for the last example in Don’s paper: 
The first line 


TRUTH_TABLE(p,q, r, (((p OR q) imp r) and not r) imp (not p and not q)) 
displays as: 

TRUTH_TABLE(p,q, r, (((p v q) -» r) a r) -» (-. p a q)) 


true 

true 

true 

true 

true 

true 

false 

true 

true 

false 

true 

true 

true 

false 

false 

true 

false 

true 

true 

true 

false 

true 

false 

true 

false 

false 

true 

true 

. false 

false 

false 

true 


Not all operators can be taken from the mathematical toolbar, some must be written 
explicitly: 

[p xor q , p iff q]= 
displays and simplifies as follows: 

[p V q, p « q] = [(-. p A q) V (p A q), (-. p v q) a (p v q)] 


Although Don’s contribution reads about a tool for the TI-89 it is no problem at all to run this 
program on the other Tl-devices, too(TI-92, Voyage 200). See V 200 screenshots: 


■' FIt Y V 

Truth Tables Wars 


truth-t <imp<a xor b.not c>> 

TRUTH FiliD RUTD FUhK 0/50 
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Usually I don’t reprint Tl-programs but I assume that many of you don’t have Graph Link for 
the TI-92/V 200 in the computers any more. Don’s program is very sophisticated from some 
points of view. He needs some auxiliary functions collaborating with two programs. The 
functions are (in alphabetical order): 

del el em ( 1 i st , i ) returns 1 i st with the ith element removed. 

del el em(l st , i ) 

Func 

augment (1 eft (1 st , i -1 ) , mid (1 st , i +1 , dim(l st) -i ) ) 

EndFunc 


i f f (a , e) and i f f (a , e) define logic equivalence and implication 

iff (a.e) 

Func 

(not a or e) and (a or not e) 

EndFunc 


imp(a.e) 

Func 

not a or e 
EndFunc 


memberg (1 i st/mat , e) returns true if e is an element of list/mat, else false 

memberq (1 i st , e) 

Func 

Local i : I f getType(l i st) = "MAT" : mat >-1 i st (1 i st)->l i st : For 
i ,1 , dim(l i st) : If when (e=l i st [i ] .true , fal se , fal se) : Return 
true : EndFor : Return fal se : EndFunc 






p 28 


MacDonald Phillips: Truth Tables on the TI-89 


DNL 98 


RmDup(l i st) returns list with duplicate elements removed - Bhuvanesh Bhatt 

rmdup(l i st) 

Func 

Local i:For i , dim(l i st) , 1 , '1 : If memberq (1 eft (1 i st , i -1 ) , 1 i st [i ] ) 
: del el em (1 i st , i )-*■! i st : End For : 1 i st : End Func 


subst ( u , o , n ) substitutes elements of n for elements in o 

subst (u , o , n) 

Func 

Local a,d,i 
d i m ( o ) ->d 
For i , 1 , d 
o[i ]=n[i ]->a 
u | a->u 
EndFor 
Return u 
EndFunc 


VarLi st (xpr) returns all variables in expression xpr - Bhuvanesh Bhatt 

varl i st (xpr) 

Func 

Local tmp 

If getType(xpr) = "MAT" :mat«-l i st (xpr)-^xpr : If getType(xpr) = "LIST" 
Then : If di m (xpr)=1 : xpr [1 ] -»xpr : End If 

If when (part (xpr) =0 , true , fal se , fal se) : Return {xpr}:If 
part (xpr) =1 : Return rmdup(varl i st (part (xpr , 1 ) ) ) : 

If when(part(xpr)>= 2 , true , fal se , fal se) 

Then : augment (varl i st (part (xpr , 1 ) ) , varl i st (seq (part (xpr , T T ) , T T , 2 , 

part(xpr) ) ) )->tmp : Return rmdup(tmp) :EndIf 

EndFunc 


menu ( ) generates a specific menu bar on the top of the screen with drop down menus 

menu ( ) 

Prgm 

Custom 

Title "Truth Tables" 

Item "Truth_T(" 

Item " not " : Item " and " 

Item " or ":Item " xor " 

Item "IMP(" : Item "IFF(" 

Title "Vars" 

Item "LogicExp" 

Item "p" : Item "q" 

Item "r" : Item "s" 

Item "t" : Item "u" 

Item "v" 

EndCustm 

CustmOn 

EndPrgm 
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truth_t(xpr) is the main program 

truth_t (xpr) 

Prgm 

Local i _ , j _ , p_ , 1 _ , m_ , k_ , v_ , n_ , o_ , q_ , c_ , md_ , xp_ 

xpr->l ogi cexp 

If when (xpr=true , true , fal se , fal se) Then 
Disp "Tautology" 

Pause 

Cl rIO : Di spHome 

Elself when (xpr=fal se , true , fal se , fal se) Then 
Disp "Self-Contradiction" 

Pause 

Cl rIO : Di spHome 
Else 

varl i st (xpr)-»-v_ 

Sort A v_ 

augment (v_, {exp} )->v1 :dim(v_)->n_ 

2 A n_->i_: { 

For j , 1 , n_ 

i_/2 A j_->p_ 

seq ( x_ , y_ , 1 , p_) ->1 _ :seq ( 1 _ , x_ , 1 , 0 , ' 1 ) ->-1 _ 
mat>-list(l_)->-l_ 

If when (j_>1 , true , fal se , fal se) Then 

For k , 1 , j_-1 

augment (1_, 1_)->1_ 

EndFor 

Endlf 

augment (m_, l_)->m_ 

EndFor 

m 1 1 =t and 0=f-»-m_: (1 i st >-mat (m_, i_) ) T -»-m_ 

m_jt=true and f=false^m_ 

{}-o_ 

For j_, 1 , i_ 

subst (xpr , v_, mat »1 i st (m_[ j_] ) )->-q_: augment (o_, { q_} )->-o_ 
EndFor 

1 i st«-mat (o_, 1 )->o_: 1 i st«-mat (vl , n_+1 )->v1 
augment(v1 ; augment (m_, o_) )->tt 
subst (tt, {true, false} , {t ,f})->-tt 
cei 1 i ng (i_/5)->c_: mod (i_, 5)->-md_ 

For j , 1 , c_ 

ClrlO 

Disp tt[1] 

If j_=c_ and md_>0 Then 

For k , 1 , md_ 

Disp tt [ ( j_-1 ) *5+k_+1 ] 

EndFor 

Else 

For k , 1 ,5 

Disp tt [ ( j_-1 ) *5+k_+1 ] 

EndFor 

Endlf 

Pause 

EndFor 

Cl rIO : Di spHome 

Endlf 

EndPrgm 
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Function varlist is of special interest because it extracts all variables in use of the given logic 
expressions. 

[1 ] http://www.technicalc.org/bbhatt/ 

This is a very rich resource for TI-89/92 & V 200 including many useful links to other 
websites 

But how is it with Tl-Nspire? TI-NspireCAS has neither built in the truth_table function nor 
can it run the TI-89 (92, V 200) program. It is necessary to write an own program. See how it 
works: 


There are two screen shots in handheld view: 




followed by screen shots from the PC: 


p 

Q 

r 

"imp (imp (p or q,r) and ] 

tine 

hue 

true 

false 

true 

hue 

false 

hue 

true 

false 

true 

false 

true 

false 

false 

true 

false 

true 

true 

false 

false 

hue 

false 

true 

false 

false 

hue 

hue 

false 

false 

false 

hue 


Done 


truth _/Q 


Enter variables: p,q,r 

Enter Boolean expression: imp (imp (p or q,r) and not r,not p and not q) 


P 

Q 

r 

"imp (imp (p or q,r) and nc 

hue 

true 

hue 

hue 

true 

hue 

false 

hue 

hue 

false 

true 

true 

hue 

false 

false 

true 

false 

hue 

hue 

true 

false 

hue 

false 

true 

false 

false 

hue 

hue 

false 

false 

false 

hue 


0.7“ 

tr 

b 






Enter Boolean expression: imp(p and s,r or q) 


OK : Cancel 


- 
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We will solve an example from [2]: 

Let p, q, rdenote the primitive statements given as: 

p: Roger studies. 

q\ Roger plays tennis. 
r Roger passes the maths exam. 

Now let pi, p 2 , P 3 denote the premises: 

p-\\ If Roger studies then, he will pass the exam. 
p 2 : If Roger doesn’t play tennis then he’ll study. 

p 3 : Roger failed the exam. 

We want to determine whether the argument (p 1 ap 2 Ap 3 j -» q is valid. 

We rewrite the premises as: 

Pi: P^r 
p 2 ; -* P 

Ps: ~^r 

Finally we generate the truth table for the implication: 

((p * r) a (-.g > p) a \P) ♦ g 


0 

Done I 

truth _t{) 

Enter variables: p,q,r 

Enter Boolean expression: imp (imp (p or q,r) and not r.not p and not q) 
p q ?' 11 imp (imp (p or q,r) and not r.not p and not q) M 


true 

true 

true 

true 

true 

true 

false 

true 

true 

false 

true 

true 

true 

false 

false 

true 

false 

true 

true 

true 

false 

true 

false 

true 

false 

false 

true 

true 

false 

false 

false 

true 


Done 


0 


The truth table presents a tautology. So we can say that (p 1 a p 2 a p 3 ) -> q is a valid 
argument. 
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KITETSMB for NspireCAS (and for DERIVE, of course) 

Michel Beaudin (Montreal, CAN) and Josef Bohm (Wurmla, AUT) 

You might remember the Integration of Piecewise Continuous Functions contributed by Michel 
Beaudin, Frederick Henri and Genevieve Savard in DNL#91. The Canadian colleagues provided a 
collection of utility functions in several KITETS... files for the TI-NspireCAS. We must not forget 
Chantal Trottier who is author of ETS specfunc.tns. All these tools can be downloaded from the 
ETSMTL website. 

Michel set an update of his papers on the web (in French, of course). One paper is a list of the func- 
tions together with their syntax provided for some of his courses at the ETSMTL. The other shows 
some examples how to apply the functions. 

It was my idea not only to combine both papers in an English version but also add examples for all 
functions appearing in the list. What you can find in the following is the first part of this list together 
with appropriate examples. Hope that you enjoy it and that you will find it useful, Josef. 

(Needless to say, that a ZJET^/FE-appendix will not be missed.) 


newton(f,x,a,n ) 


n times application of Newton’s method for solving an 


► *Kit_ETS_MB ■ 


1.1 | 1.2 | 13 
©VERSION DU 20 AVRIL 2015. 


hit_ets_mb ^newton fx-cos (xj.x, 0.7,5) 


mm 


S 


equation of form j(x) = 0 with a as initial value (first es- 
timation). 


0.7 

0.739436 

0.739085 

0.739085 

0.739085 

0.739085 




fixed _point(fx,a,n ) 

n times application of the fix point method for solving an 
equation of form x = f(x) with a as initial value (first es- 
timation). 

I transposed the result in order to receive a column which 
is easier to be read. 

We observe that this method converges much slower than 
Newton’s method. I present the result after 20 iterations. 
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leftsum(f,x, a, b, n ) 
rightsum(f,x, a,b,n) 
midsum(f,x, a, b, n ) 

give the respective sums using n subdivisions for estimat- 

b 

ing the integral J / (x)dx. 


Unfortunately NspireCAS CX is unable to find a closed 
form for some sums (eg. for the trig functions). 

We can calculate the sums for large values of n numeri- 
cally. 


H 1.1 1 1.2 1 1.3 wh 

ETS_MB ^ <1Q| 

& solve (x-cos (x) = Q,x) 

x=0. 739085 

- 

Ieflsttm[z^ t z, 0 , 1/7 ) 

2 n 2 -Z n+1 



6-T7 2 


/ 2 1 
2-77^-3-77+1 

1 


"-”1 6 n 2 1 

3 

1 


■ 

s | 


i.i 


X 2 I 

lefisnm (sin (. z),z , n, 2 ■ n.n ) 


3 0 


n - 1 


7 = 0 


f , j 7- n \ 
sin 1 1 jc 

t l " 


0 



IB 1 1 


► *Kit_ETS_MB^ 


Nspire produces the sum and returns it in exact form by 
presenting all summands. 

The closed form for the left sum is given by 


71 

— cot 
n 


71 


v2 nj 


Unlike NspireCAS is DERIVE able to find the closed 
forms for the approximating sums of trig functions as you 
can see below. 

[F(x) := 5IN(x), a := tt ? b := 2 -tt] 

[LSUM_VAL(n) , RSU^LVAL(n), MSUM_VAL(n)] 


r(n ): =ngfttsum (sin (z) d z d n d 2 ■ n d n) 
m(n)\ =midsum (sin [z),z, n, 2 ■ n, n ) 
/( 2 0 □ ) 


Done 

Done 

Done 


/ 39- 


19- n \ 


' 37 3 

\ 2O0 ) 


, ioo 1 


, 200* 

100 


100 


100 


- 


""fe) "“{it) 




lim [LSUM_VAL(n), RSUM_VAL(n) , MSUM_VAL(n)] 
n-«i 


[- 2 , - 2 , - 2 ] 



1 1 2 j 1 3 u 

*Kit_ETS_MB ^ 

ma 

r 

[ 200 1 

l 1O0 j 

\ 200^ 

3 


100 

100 

100 


7 ( 200 ) 


-1.99996 


r 

(300) 


-1.99998 


777 ( 500 ) 


-2. 



2- ji 

sin(z) d z 


-2 

1 

- 

n 


L_ III 
>|| 


trap(f,x,a,b,n) 

Trapezium method for approximating a defined integral. 


1.3 


y 11 l__l 

trap (x " J +x,x, 1 . 4 / 7 ) 


RADljlQ 


15- 19-t? z +9 


4- 77 ^ 


lim 


15- 19-77^+9 


4- 77 


285 

4 


,x J +xJd* 


285 


- 
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simpson (f,x, a, b, n) 

Simpson method for numerical integration ( n must be an 
even number). 



I don’t receive an answer in closed form even when I 
enter In subdivisions (even). 


O' 


1.1 1.2 1.3 


53.59S2 


e 4_l V + 4 e 10 + lj 


30- \e J -1/ 


n 1 1 j 1 


trap Lx +x,x, 1.4/j, 1 


15- 19- r+9 


2, n = 


4- n 


15- (l9-rc 2 +9^ 

■ J 



2S5 

4 




I substitute for n = 20 and receive a pretty 
expression which an be approximated 


simpson ( e x ,x, 0, 4, 2 0 ) 

19 IS 17 16 

e 4 +4- e 5 +2- e 5 +4- e 5 +2- e 5 +4-e* 

B 


Entering n = 20 directly in the function gives the sum 
explicitly with all its summands. 



Approximation results in the same value. 
Let’s compare with DERIVE : 


x 

SIMPS0N(e p x f 0, 4 P 2-n) 


4 2/n 1/n 

(e - !)■ (e + 4-e + 1) 


2/n 

3 -n- (e - 1) 


x 4 

lim SIMPS0N(e , x, 0 P 4, 2-n) = e - 1 

n-wo 
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ratio _test(t,n) 


Gives the lim- 


tin) 


for testing the convergence of the 


t(n + 1) 

oo 

sum of the series ^\(«). The series is convergent if the 

n = 1 

limit is less 1 , it is divergent if greater 1 . If the ratio test 
result is equal 1 then we cannot decide (and possibly try 
other tests). 


When I read about the ratio test I remembered that I had 
learned about a root test, too - some decades ago. 

The sum of a series converges if lim^ < 1. 

«-> 00 

Defining the function root_test{t,n ) was an easy job. 

First series is convergent, the second one is not. 


Ji n’i i2l i~TT| 


n\ 

ratio _test\ — ,n 

W" 


ratio _test\ 


x n 2 n ~ l 
2 n-1 


-,n 


PADflJQ 


-1 


2\x\ 


ratio Jest(pCv(2 n,n ) x n ,n t 


4 \x\ 


n ! . 

2 ^ — ~ is convergent. 

n = 1 n 


The second series is convergent for 
|x| < 14 and the third one is convergent for 

M < ‘A 


n 1 1 1 1 1 12I 


1.131 


*Kit_ETS_MB ^7 


Define Lib Pub root_test{t,n)= lim \ \[tj 
n 

^ Fertig 


root_test\\l — ,n 

n I 


root testy 


-1 


1 I 

2 


\2 ■ n\ 


ang_vect(u ,\ ) 

returns the angle (between 0 and n) between the two vec- 
tors u and v(u,ve R 2 or R 3 ). 


1 EKDEKBElEh *k*-ets_mb^ 

pad Bill 

ang_2_vect([ 2 3],[-4 l]) 


3 

• i 5 

•[221 \ ji 


sin' 1 — 

+ 


* \ 

221 / 2 


ang_2_vect([ 2 3 ], [ - 4 l]) 

1.91382 


ang_2_vect([2 0 o]^[ 0 0 -3]) 

n 

2 


ang_2_vect([ 2 2 o],[l 0 - 3 ]) 

1.34528 

■ 


eI| 


atan(y,x) 

returns the angle 9 (-n < 9 < n) of the vector [x,y]. 
Notice the order of entering the coordinates. 


Hi 11 1 1 1 2 1 1 1 3| 


► *Kit_ETS_MB r 


DEG <SO 


I atari 2(2, -3) 


2.55359 

- 

© change to degree mode 



atan2(y, l) 


45 


atan2^2 t -^) 


180-tan M |— -j 


a tan 2{2 t -2 > ) 


146.31 

1 



*Kit_ETS_MB ^ 


RADiJJO 


a tan 2(l , l) 
atan2^2,-2) 

0 ^ 2 ( 2 , - 3 ) 

© change to degree mode 


jr-tan‘|— 
2.55359 


In radiants (above) and in degree mode 
(left). 
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ptcri(f ,\ ) 

returns the matrix containing the critical points of func- 
tion / of two variables v = [x,y].f should be a polynomial. 

Given is f{x,y) = 10 x 2 y - 5x 2 - 4 y 2 - x 4 - 2 y 4 . 

The critical points can be found by solving the system of 
equations V/ (x, y) = [f x , f y ] = [0, 0] . 

nature(f,\,\ o) 

The nature of a critical point vo = \a,b] can be decided by 
evaluating the determinant of the Hessian matrix D to- 
gether with the value of the 2 nd partial derivative f xx in the 
respective point vo. 

This function returns [D(a,b),f xx (a,b)\. 

In the example we would like to find the nature of the 
fifth point in the list of critical points. 
As D > 0 and f xx < 0 / has a relative maximum at 
(-2.644,1.898) with function value 8.496. 

The 3D-plot of the function confirms the calculations. 

It is in fact a global maximum. The maximum point is 
realized in the graph as a sphere with its centre at (-2.64, 
1.90, 8,50) and radius 0.3. 


□ l 11] 1 12|l 13j[J 

_ETS_MB ^7 

de& <BO 

Kxj^): = 10x 2 X“5x 2 

! -4x 2 -x 4 - 

•2/ 


* 



Done 


ptcn{f(y,y),[x y])^m 






2.64422 

1.89838 




0.856657 

0.646772 




0 

0 




-0.856657 

0.646772 




. -2.64422 

1.89838 . 


, 

1 



■|| 


Din 1 1 2 1 1 t IXS 

_ETS_MB ^7 

DEsflUEJ 


0.856657 

0.646772 


3 


0 

0 




-0.856657 

0.646772 




. -2.64422 

1.89838 



nature(/(icy),[x y],m[s]) 





[2488 72 

-55.9354] 


Am[ 5,1 W5,2]) 


8.49586 


r -. 

* 



1 




M 


I hst 


113 114 I 


*Kit_ETS_MB ^ 


DEGa CHO 



lagrangel (f,g, v ,A) 

simplifies to a list of equations which can be used to ap- 
ply the method of Lagrange multipliers with one con- 
straint. / is the objective function, g the constraint (in 
form g = 0), v is the vector of the variables and k the mul- 
tiplier. 

Example: 

What is the rectangle of maximum area that can be in- 
scribed in the ellipse given by 16x 2 + 9 y 2 = 144 ? 

The vertex of the rectangle where x, y > 0 is at 


*Kit_ETS_MB ^ 


fll.11 1 1 . 1 2 | l 1 . 1 3] 

© First example for Lagrange multipliers 

la grange 1 (4 x ■ y, 1 6 • x 2 + 9 • y 2 - 1 44, [x y], a) 
{±y-22A.x,±x-18Xy,16x 2 +9y 2 -li> 
© Solve the system 

zeros({ 4 y-32- Ax,4 x-18- A y, 16- x 2 +9- jp 

lL. 2-ji - 


s 


3V2 


,2\fl 


and the maximum area is 24. 


We can see the complete calculation on the PC-screen. 


1.11 1.12 




*Kit_ETS_MB ^ 


-3- / 2 


2 {2 -i 

2 6 

- 2 ji a 
2 6 


3 ■ J 2 

4 xy\x= and y=2- \f2 


24 
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Example: 

The Cobb-Douglas production function for a special product is given by / (x,y) = 100x 3/4 y 14 where x 
represents the units of labor (at € 150 per unit) and y represents the units of capital (at € 250 per unit). 
The total expenses for labor and capital is limited to € 50.000. What is the maximum production level 
for this good? 



50 units of capital and 250 units of labor guarantee the maximum production of 16.718 product units. 
Economists call 0.334 the marginal productivity of money (i.e. each additional spent € results in 0.334 
additional product units. 
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Iagrange2(fg, v,X,jli) 

simplifies to a list of equations which can be used to ap- 
ply the method of Lagrange multipliers with two con- 
straints and three variables. / is the objective function, g 
the vector of the constraints (in form g = 0), v is the vec- 
tor of the variables and X and ju are the multipliers. 



temp (xj/,z) =20+2- x+2y+z 2 Done 

g =x 2 +y 2 +z 2 -ll:h: =x+y+z- 3 *+X+z- 3 

la grange 2(temp (x^z), [g /;],[x y z\xp) 

{ -2 A x-p+2, -2 ■ A y-p+2, -2 • (a- l) • z-p,x 2 J 
ZQTOs{lagrange2{temp(^ J y t z) J [g a],[x y z\ v 


3 

-l 

1 0 

~ nr* 

~ nr.* 

l rr *\ - 


Example: 

temp(x,y,z) = 20 + 2x + 2y + z 2 describes the temperature at each point on the sphere with its centre in 
the origin and radius Vl 1. What are the extreme temperatures on the intersection curve of the sphere 
and the plane x+y + z= 3? 



The maximum temperature is 30.33 and the minimum temperature is 25. 

(Note: The examples with the Lagrange multipliers are from Larson, Hostetler, Edwards : Multivariable Calcu- 
lus.) 
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curl{ F,v) 

returns the rotation (in German rctf(F,v)) of a vector field 
of three components. 

Find the curl of F = [x 2 z, -3xz 2 , xyz ] at (2,-1, 4) 

Find the rotation of [x sin x, y sin y, z sin z]. 

(Vector fields with curl = zero vector are called irrota- 
tional. F is conservative iff its curl is 0.) 


*Kit_ETS_MB t 


!»*» 10 O 


curl^x 2 z -3 xz 2 X7’zl[x y z]) 

\_ 7 xz x 2 -yz -3 -z 2 ] 

7 xz x 2 -yz -3 z^]|x=2 andx=’l andz=^ 

8 -48] 

y z\) 

0 0 o] 


[56 8 -48 
ci/r/([x sin(x) y sin(y) z- sin(z)],[x y z]) 


potential^ F,v) 

returns the potential of a vector field. 

Check always the answer by calculating the gradient of 
the result. 


DTTTT 


First check if the curl of F = 
find its potential. 


1 -x ! 

? 2 1 

y y 


= 0 and if so. 


1.12 


curl\ 


Fm 


*Kit_ETS_MB ^ 


1 -x 

- — 2z-l 


I* y z] 


PADiJJQ 


[o 0 o] 


potential 


ill 

-x 

r i\ 

- 

— 2- z-l 

.1* y z\\ 

P 

y 2 

1 


(z-l) 

gradip (z-l),[x y z]) [o 0 2 z-lJ 


0 


PAD i0Q 


Pl-1lll.12l l.13 

gradi?- (z-l),[x y z]) [o 0 2-z-lJ 

© This is not the complete potential . It is 
— +z (z - 1) +constant 


grad\-+z- (z-l)+c,[x y z 

\y 


i - 2- z-l 
' y 2 


a 


Find the potential of G(x, y, z)=[y 2 z 3 , 2 xyz 3 , 3xy 2 z 2 ]. 

It is easy to see - without CAS - that it is given by xyz . 

Verify that the gradient of the result gives G from above. 

It is important that certain line integrals can be calculated 
by means of the potential. 

Calculate the line integral in this vector field along the 
curve given by r(t) = [5 1, 3 1 - 1, t\ 1 < t< 4. 

Then apply function int_curv contained the library. 


f 


1 11 1 11 2 1 1. 1 3 

2 


*Kit_ETS_MB ^ 


7 7 9 9 1 

z 2 x y z Z x y z ]-+g 


[ y 2 z 3 2-X7 Z 3 3-X7 2 -z 2 J 


[x y z]-v 
curl(g t v ) 
potential(g,v ) 

xy 2 z 2 ->pot{x,y t z) 


Lx y z\ 

[o 0 o] 

2 3 

xy z 

Done 
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int_curve(F,x,r,t,a,b) calculates the line integral in the 
vector field F (variables v) along the parameter curve r(f) 
with a<t<b. 


The vector field H(x, y, z) = — -[-y,x,0] is not con- 

x I 2 +y 2 

servative because its line integral over a certain curve 
does not pass the origin. But its rotation (curl) is the zero 
vector 0 . 

The domain of H are the points in space except those on 
the z-axis. We choose as closed path the circle x 2 +y 2 = 1, 
z — 0 (which is a circle in the xy-plane). 

Show that the line integral is 2n. 

We find the integral without applying int_curv. 


*Kit_ETS_MB ' 


rad <50 


& curl(h(x A y,z) t v) 
[cos(/) sin(/) o] ->r{t) 


[o 0 o] 

Done 


2 - tc 


dotP[/7(cos(/),sin(/), o), ^ (/'(/)) ) d/ 


2 71 


0 



I add three more examples: 

(1) Integrate the vector field [y 2 , 2 xy - along [cos(f), sin(/)], 0 < t < n/2. 



This was an example in the xy-plane. Let’s do a similar problem in space: 

x v 1 

(2) Find the work done by the force field F(x,y,z) = — 5 — on a particle as it moves along a 

4 2 4J 

helix like space curve given by r(f) = cos(f), 2sin(^), — , 0 <t< An. 
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( 3 ) 


The components of a force G are given by 


x + y-3z a-x + b-y + c-z d-x + e-y + f-z 
(x 4- y 4- z)^ (x 4- y + (x 4- y 4- z)^ 


Find the coefficients a, b, c, d, e and / such that the force has potential u and find G and u. 
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The last example cannot be solved by means of Michel’s tools. We integrate fx wrt x (because it is the 
partial derivative of the possible potential. We integrate fy and fz wrt to y and z and perform a com- 
parison of the coefficients. 


newtons(u,x,vO,n) applies Newton’s method to find a 
solution of the system u of more variables (given in vec- 
tor v) with initial solution vector vO performing n itera- 
tions. 

The calculation process can last a longer time. 

It can also happen that one is not successful with the cho- 
sen initial vector. Then try another one. 

Example: 

Solve the system 


2 2 

x + y ■ 


■z 2 =-5 


2y + z 2 = 20 
3x + z = 40 


► *Kit_ETS MB c 


RAD <SO 


[-4 47 18 61 41.6S] 


f +z 2 - 20 3x 2 +z-4o].[x y z^[o 0 0I5) 
"Error: Argument must be a variable name." 

*+z 2 -20 3 x 2 +z-4o],[x y z}[ 0 0 o],l)' 


0 0 0 . 
undef undef 00 


newton J^x^ +y* -z^ +S 2 y+z^ -20 3x 2 4* 0 


► *Kit_ETS_MB^ 


0. 0. 0. 
undef undef 00 


Vz 2 -20 3x 2 +z-4o],[x y z\[l 1 l],5, 


1 . 1 . 1 . 
5.65385 1.42308 9.07692 
3.8435 1.1411 5.51444 

3.46622 1.03447 4.38305 
3.45255 1.01998 4.24033 
3.45263 1.01973 4.23799 


0 


lim_mat( mat,v,v0) returns the limit of mat when v tends 
to vO. mat can be a scalar vector field, a vector or a ma- 
trix. 


1.6 1.7 


► *Kit_ETS_MB c 


P 15 1 

lim_ma&Z x y x 2 >’ 2 -x],[x ^],[l 2]) 

[6 1 3] 

lUM) 

[(a+b) 2 Z (a+b) 2-ff+Z?] 


/7/77_/77£7/([x 2 3‘X X+Ct 


1im_rnat 


xy x 

\y 2 xy] 



6 4 
9 6 


M 


newtons needs the ITERATES-command which is not a Nspire built-in function. It was not so easy to 
transfer this great DERIVE - function to the Nspire-technology. In early TI-92 days this was done by 
Terence Etchells - as far as I do remember. You may find the discussion between Michel and me in- 
teresting on this issue which is published in D£7?/Fi?-Newsletter #90. 


The next screen shot shows the PC-screen containing some more examples of applying newtons and 
lim mat as well. 
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newtons I 


(sin (x) ) 2 2 ■ x+>' 2 -4 , [x y\ [ 1 1 ], 5 




([(smfx)) 2 -^ 2x+y 2 -i\[x y\[ 


1.50439 0.995597 


newtons ([x 3 - 2 ■ x 2 +y - 1 0 l-x J +2z 100+x^-y 2 +z 2 -10G],[x y z],[ 1 1 l], l) 


newtons [\x 3 - 2 - x 2 +y - 1 0 l-x J +2z 100+x^-> ,2 +z 2 -100 


\i x y 4 L 1 1 43) 


1. 1. 

1.41477 1.08523 

1.50108 1.00235 

1.5044 0.995609 

1.50439 0.995597 
1.50439 0.995597 

[-3.21936E-7 -0.000007] 

1 . 1 . 1 . 

^5.2 15.2 6.3 

1 . 1 . 1 . 

5.2 15.2 6.3 

3.77992 9.13151 12.2054 
2.92322 8.33689 8.14282 


4 


hm_mat\\x' J -2 - x 2 +>'-10 l-x 3 +2z 100+x^-> ,ji, +z^-10o],[x y z\\_2. 


923 8.337 8.143]) 

[6.22305 -7.6879 69.8016] 


newtons 


([ x 3 -2 x 2 +>'-10 1 - x 3 + 2 z 100+x| 4 ->' 2 +z 2 -10o],[x y z}[3 8 S^s) 


3. 8. 8. 

2.51596 8.26064 6.46543 
2.36098 8.11734 5.99154 
2.34531 8.10186 5.94932 
2.34516 8.10169 5.94893 
2.34516 8.10169 5.94893 


J 0 


Michel sent a rich collection of more utility functions and examples how to apply them. We will pub- 
lish his materials in the next newsletters. 

He also sent a DERIVE file containing some of his functions. I accomplished this collection of tools 
and tried to use the same function names and parameters as well. You can find the respective file on 
the next two pages. 

I added three examples from above performed with DERIVE. 


References: 

[1] Gunter & Kusmin , Aufgabensammlung zur Hoheren Mathematik, VEB Deutscher Verlag der 
Wissenschaften, Berlin 1964 

[2] T. Arens a.o ., Mathematik, Spektrum 2009 

[3] Larson, Hostetler, Edwards, Multivariable Calculus, Houghton Mifflin Company 2002 
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Michel’s Tools prepared for DERIVE 


ki tETS_DERIVE . mth 


NEWTON(f , x, xO, n) = ITERATES 


— f a x a xO , n 
dx 


fixed_point(f a x a xO , n) := ITERATES (f, x a xO „ n) 

n b - a 

RIGHTSUMfg, x a a, b a n) := Z ( lim g)- 

i=l x-*a + i ■ (b - a)/n n 


n - 1 b - a 

LEFTSUMCg, x, a, b, n) := Z ( lim g) 

i=0 x-*a + i - (b - a)/n n 


n b - a 

MIDSUMCg, x, a, b, n) := Z ( lim g) 

i=l x-»a + (i - l/Z)-(b - a)/n n 


TRAPfg, x, a, b , n) := 


(b - a)- 


(( 


n 

Z lim 

i=l x-n - (b - a)/n 



n - 1 

+ Z lim 

i=0 x->i ■ (b - a)/n 


g 

-i- a J 


2-n 


SIMPSONCg, x , a, b, n) := 

If MQD(n a Z) t 0 

"n pair S.V.P. ! ! ! " 

(Z ■ MIDSUMCg , x, a, b a n/Z) -f TRAPfg, x, a. b, n/2))/3 


ratio_test is implemented in DERIVE 

1/n 

root_test(u a n) := lim u 
n-K» 


ang_Z_vect(y , v) := ACOS 


u ■ v 1 

M-M J 

atan2(y,x) is implemented in DERIVE 


ptcriff, v) := SOLUTIONS (3 (f , v ) a 3(f, v ) , |"v , v 1) 

1 Z L 1 zj 


Z 

hessian (f s v) := 3(f, v a 2)-3(f a v a 2) - 3(3(f a v ), v ) 

1 Z 1 Z 


nature(f, v, vO) := Him hessian(f, v) s lim 3(f, v a 2)1 

L v^vlJ v^vO 1 J 

NEWTONS(u a v, vO ,r) is implemented 

lagrangelff, g, v := [x, y] a X) := APPEND(GRAD(f - X-g, v) s [g]> 

lagrangeZCf, g, v := [x, y, z] a X, p) := APPEND(GRAD(f - X-g - p-g a v), 

1 Z 


g) 
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GRAD(f,v) a DIV(F,v) a LAPLACIAN(F , v) and QJRl(F,v} are implemented, see Online Help. 
jacobian(f, v) := VECTOR (VECTOR (3 (f a v ) 3 i s DIMCv)), ] a DIM(f)) 

i ^ 

lim_mat is not necessary, LIM-function is sufficient for all cases! 

1 


potential (f, v) 


:= T ( lim 

J v-^t_ ■ > 


f)-v dt_ 


0 


u 

i 


int_curv(ch, ve , co , pa, a, b) := | ( lim ch)- 

ve-^co d pa 
a 

Three examples from above worked with DERIVE : 


co dpa 


Z Z Z 4 4 

ptcriflOx - y - 5 - x - 4-y - x -Z-y a [x, y]}= 


L 

— L 


0.3566568714 0 .64677Z1990 

-0 .3566563714 0 . 64677Z1990 

Z . 644ZZ4301 1 . 893334431 

-Z . 644ZZ4301 1 . 398334431 

3 . 90Z0ZZ956 - 1 -Z . 545156630 

L - 3.90Z0ZZ956 -l -Z. 545156630 J 


z z z z 

_temp(x :s y a z):=Z0 + Z-x + Z-y + z a g:=x + y + z -11, h :: x + y 


+ +z - 3] 


SGLUTI0NS(la.grange2(temp(x, y, z), [g, h], [x, y, z]. A, p), [x, y, z, A, p]> 
-1 3 10 2 

3 -1 102 

2 2 8-/3 

3 3 9 


2-/3 2-V3 4-/3 

+ 1 + 1 1 


2-/3 2-/3 4 - J3 2 8-J3 2 

1 1 + 1 + — 

3 3 3 3 9 3 


r(t) := 


int_curv 


int_curv 


COSCt), SIN(Z-t), 

x y 1 

4 Z 4 . 

x y 1 

4 Z 4 . 


[x, y, z], r(t), t, 0, 4 -jt 


[x, y, z], r(t), t, 0, 4 -jt 


JT 

z 


